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Abstract 

The decays B^!^ — >• J/i|j Dg^ and B+ — )■ J/4) D*+ are observed for the first time using 
a dataset, corresponding to an integrated luminosity of 3fb~^, collected by the 
LHCb experiment in proton-proton collisions at centre-of-mass energies of y/s = 
7 and 8 TeV. The statistical significance for both signals is in excess of 9 standard 
deviations. The following ratios of branching fractions are measured to be 



g(B+^J/r|;D:+ ) 



2.90 ±0.57 ±0.24, 
2.37 ±0.56 ±0.10, 



where the first uncertainties are statistical and the second systematic. 
The mass of the B+ meson is measured to be 

= 6276.28 ± 1.44 (stat) ± 0.36 (syst) MeV/c^ 

using the B+ — J/\|j D+ decay mode. 
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1 Introduction 

The B+ meson, the ground state of the be system, is unique, being the only weakly decaying 
heavy quarkonium system. Its lifetime [T|[2] is almost three times smaller than that of 
other beauty mesons, pointing to the important role of the charm quark in weak B+ decays. 
The meson was first observed through its semileptonic decay — t- J/il;£"'"V£X [3]. Only 
three hadronic modes have been observed so far: B+ — t- J/i[)7r+ |4|, B+ — t- J/ip Tr+yr+Tr" ^ 
and B+ ^•il;(2S)7r+ |6]. 

The first observations of the decays B+ J/i[)D+ and B+ — > J/vIjD*"'" are reported in 
this paper. The leading Feynman diagrams of these decays are shown in Fig. [T| The decay 
B+ — )■ J/i[)D+ is expected to proceed mainly through spectator and colour-suppressed 
spectator diagrams. In contrast to decays of other beauty hadrons, the weak annihilation 
topology is not suppressed and can contribute significantly to the decay amplitude. 




Figure 1: Feynman diagrams for B+ — )• J/\);D+ decays: (a) spectator, (b) colour-suppressed 
spectator and (c) annihilation topology. 

Assuming that the spectator diagram dominates and that factorization holds, the 
following approximations can be established 

^ r(B+^J/x|;D+) ^ F (B ^ D*D+) 
/^D+/.+ - r ^ ^ F (B ^ D*7t+) ' ^ ' 

^ F(B+^J/x|;D:+) ^ F(B^D*D:+) 

- F(B+^J/i|;D+) F(B^D*D+) ' ^ ' 

where B stands for B"^ or B and D* denotes D or D*~. Phase space corrections amount 



to (9(0.5%) for Eq. (la) and can be as large as 28% for Eq. (lb), depending on the relative 
orbital momentum. The relative branching ratios estimated in this way, together with 
more detailed theoretical calculations, are listed in Table [TJ where the branching fractions 
for the B — D*D+ and B — )■ D*7r+ decays are taken from Ref. [l]. 

The analysis presented here is based on a data sample, corresponding to an integrated 
luminosity of 3fb~ , collected with the LHCb detector during 2011 and 2012 in pp collisions 
at centre-of-mass energies of 7 and 8 TeV, respectively. The decay B+ J/i|; 7t+ is used as 
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Table 1: Predictions for the ratios of meson branching fractions. In the case of 7^j-,*+^j-,+ 
the second uncertainty is related to the unknown relative orbital momentum. 





2.90 ±0.42 
1.58 ±0.34 



1.3 
2.6 
2.2 
1.2 



2.20 ±0.35 ±0.62 
2.07 ±0.52 ±0.52 

3.9 

1.7 



Eqs. Q with B° 
Eqs. Q with B+ 

Ref. g 

Ref. 8 

Ref. 9 



Ref. 10 



a normalization channel for the measurement of the branching fraction i3(B+ — )■ J/i|;D+). 
In addition, the low energy release (Q- value) in the — )■ J/i[) mode allows a determi- 
nation of the mass with small systematic uncertainty. 



The LHCb detector [TTJ is a single-arm forward spectrometer covering the pseudorapidity 

range 2 < rj < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes placed downstream. The combined tracking system has 
momentum resolution Ap/p that varies from 0.4% at 5 GeV/c to 0.6% at 100 GeV/c, and 
impact parameter resolution of 20 |J.m for tracks with high transverse momentum. Charged 
hadrons are identified using two ring-imaging Cherenkov detectors. Photon, electron and 
hadron candidates are identified by a calorimeter system consisting of scintillating-pad and 
preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are 
identified by a system composed of alternating layers of iron and multiwire proportional 
chambers. The trigger [l2] consists of a hardware stage, based on information from 
the calorimeter and muon systems, followed by a software stage which applies a full event 
reconstruction. 

This analysis uses events collected by triggers that select the decay products of 
the dimuon decay of the J/i|» meson with high efficiency. At the hardware stage either 
one or two identified muon candidates are required. In the case of single muon triggers 
the transverse momentum, p^, of the candidate is required to be larger than 1.5 GeV/c. 
For dimuon candidates a requirement on the product of the of the muon candidates is 
applied, y/PTiPT2 > l-3GeV/c. At the subsequent software trigger stage, two muons with 
invariant mass in the interval 2.97 < < 3.21 GeV/c^ and consistent with originating 

from a common vertex are required. 

The detector acceptance and response are estimated with simulated data. Proton- 



2 LHCb detector 




2 



proton collisions are generated using Pythia 6.4 
in Ref, 
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13j with the configuration described 
Particle decays are then simulated by EvtGen 15 in which final state 
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radiation is generated using Photos 

the detector and its response are implemented using the Geant4 toolkit 
in Ref. 



The interaction of the generated particles with 

as described 
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3 Event selection 

Track quality of charged particles is ensured by requiring that the per degree of freedom, 
Xtr/iidf, is less than 4. Further suppression of fake tracks created by the reconstruction 
is achieved by a neural network trained to discriminate between these and real particles, 
using a requirement on the output of this neural network, Pfake < 0.5. Duplicate particles 
created by the reconstruction are suppressed by requiring the symmetrized KuUback- 
Leibler divergence fl9], Agjf, calculated with respect to all particles in the event, to be in 
excess of 5000. In addition, the transverse momentum is required to be greater than 550 
(250) MeV/c for each muon (hadron) candidate. 

Well identified muons are selected by requiring that the difference in logarithms of 
the likelihood of the muon hypothesis, as provided by the muon system, with respect to 
the pion hypothesis, A^/"ln£ [20], is greater than zero. Good quality particle identification 
by the ring-imaging Cherenkov detectors is ensured by requiring the momentum of 
the hadron candidates, p, to be between 3.2GeV/c and lOOGeV/c, and the pseudorapidity 
to be in the range 2 < rj < 5. To select well-identified kaons (pions) the corresponding 
difference in logarithms of the likelihood of the kaon and pion hypotheses [2l] is required to 
be A^/"ln£ > 2(< 0). These criteria are chosen to be tight enough to reduce significantly 
the background due to misidentification, whilst ensuring good agreement between data 
and simulation. 

To ensure that the hadrons used in the analysis are inconsistent with being directly 
produced in a pp interaction vertex, the impact parameter x^j defined as the difference 
between the of the reconstructed pp collision vertex formed with and without the 
considered track, is required to be xip > 9. When more than one vertex is reconstructed, 
that with the smallest value of xjp is chosen. 



As in Refs. 22 24 the selection of J/v|j — )■ \i^\L~ candidates proceeds from pairs of 
oppositely-charged muons forming a common vertex. The quality of the vertex is ensured 
by requiring the of the vertex fit, xlx^ to be less than 30. The vertex is forced to be 
well separated from the reconstructed pp interaction vertex by requiring the decay length 
significance, ^flight > 3. Finally, the mass of the dimuon combination is required to be 
within ±45MeV/c^ of the known J/ij; mass [I], which corresponds to a ±3.5cr window, 
where a is the measured J/t|) mass resolution. 

Candidate mesons are reconstructed in the D+ — )■ (K^K^)^n^ mode using criteria 
similar to those in Ref. [25|. A good vertex quality is ensured by requiring xlx < 25. 
The mass of the kaon pair is required to be consistent with the decay (() — i- K+K^, 
I^K+K- ~ "^4)1 < 20MeV/c^. Finally, the mass of the candidate is required to be within 
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±20MeV/c^ of the known D+ mass [1| and its transverse momentum to be > 1 GeV/c. 

Candidate mesons are formed from J/4) pairs with transverse momentum in excess 
of 1 GeV/ c. The candidates should be consistent with being produced in a pp interaction 
vertex by requiring Xip < 9 with respect to reconstructed pp colUsion vertices. A kinematic 
fit is appUed to the candidates |l26] . To improve the mass and Ufetime resolution, in this 
fit, a constraint on the pointing of the candidate to the primary vertex is applied together 
with mass constraints on the intermediate J/il) and states. The value of the J/4> mass 
is taken from Ref. jlj. For the meson the value of mp+ = 1968.31 ± 0.20MeV/c^ 
is used, that is the average of the values given in Ref. 27 and Ref. [l]. The P^r 



degree of freedom of this fit, Xat/^^df, is required to be less than 5. The decay time of the 
D+ candidate, cr (D+), determined by this fit, is required to satisfy cr > 75 |J.m and the 
corresponding signed significance, Scr, is required to be in excess of 3. Finahy, the decay 
time of the candidate, cr (B+), is required to be between 75 \mi and 1 mm. The upper 
edge, in excess of 7 lifetimes of meson, is introduced to remove badly recontructed 
candidates. 



4 Observation of B+ J/i|) D+ 

The mass distribution of the selected B^ — )■ J/x|; candidates is shown in Fig. |2] The peak 
close to the known mass of the B+ meson [l||28] with a width compatible with the ex- 
pected mass resolution is interpreted as being due to the B+ — )• J/i|)D+ decay. The 
wide structure between 5.9 and 6.2 GeV/ is attributed to the decay B+ — )■ J/il^D*"*", 
followed by D*"*" — )■ D+y or D*"*" — )■ D+tt" decays, where the neutral particles are not 
detected. The process B+ — )■ J/tJ^D*"*" being the decay of a pseudoscalar particle into 

two vector particles is described by three helicity amplitudes: A++, Aqq and A , where 

indices correspond to the heli cities of the J/i|) and D*"^ mesons. Simulation studies 

show that the J/i|)D^ mass distributions are the same for the A++ and A amplitudes. 

Thus, the J/4'D+ mass spectrum is described by a model consisting of the following 
components: an exponential shape to describe the combinatorial background, a Gaussian 
shape to describe the B+ — )• J/t|)D+ signal and two helicity components to describe the 
B+ — > J/i[>D*+ contributions corresponding to the A±± and ^oo amplitudes. The shape of 
these components is determined using the simulation where the branching fractions for 
D*+ — )■ D+7 and D*+ — )■ D+7r° decays are taken from Ref. [l]. 

To estimate the signal yields, an extended unbinned maximum likelihood fit to the 
mass distribution is performed. The correctness of the fit procedure together with the 
reliability of the estimated uncertainties has been extensively checked using simulation. 
The fit has seven free parameters: the mass of the B^ meson, m^+, the signal resolution, 
(Tb+, the relative amount of the A±± helicity amplitudes of total B+ — )■ J/i[)D*+ decay 
rate, f±-i-, the slope parameter of the exponential background and the yields of the two 
signal components, -^B+^j/\pD+ -^b+^.j/\|)D*+' background. The values of 

the signal parameters obtained from the fit are summarized in Table [2} The fit result is 
also shown in Fig. [2] 
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Figure 2: Mass distributions for selected J/i1jD+ pairs. The solid curve represents the result 
of a fit to the model described in the text. The contribution from the — )• J/\[)D*+ decay is 
shown with thin green dotted and thin yellow dash-dotted lines for the A±± and ^oo amplitudes, 
respectively. The insert shows a zoom of the mass region. 



To check the result, the fit has been performed with different models for the signal: 

and a modified Novosibirsk function [31] . 



a double-sided Crystal Ball function 29 , 30 



For these tests the tail and asymmetry parameters are fixed using the simulation values, 
while the parameters representing the peak position and resolution are left free to vary. 
As alternative models for the background, the product of an exponential function and 
a fourth-order polynomial function are used. The fit parameters obtained are stable with 
respect to the choice of the fit model and the fit range interval. 

The statistical significance for the B+ J/ip signal is estimated from the change in 

the likelihood function = \j2hi^^^, where is the likelihood of a background-only 
hypothesis and Cb+s is the likelihood of a background-plus-signal hypothesis. The sig- 
nificance has been estimated separately for the B^f — t- J/i|> D+ and B+ — )■ J/i[) D*+ signals. 
To exclude the look-elsewhere effect p2l, the mass and resolution of the peak are fixed to 
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Table 2: Signal parameters of the unbinned extended maximum likelihood fit to the J/ipD^ mass 
distribution. 



Parameter Value 



[MeV/c^] 6276.28 ±1.44 

[MeV/c^j 7.0 ± 1.1 

<+^j/^D+ 28.9 ± 5.6 

B+^j/^D*+ 2.37 ±0.56 



52 ±20 



the values obtained with the simulation. The minimal significance found varying the fit 
model as described above is taken as the signal significance. The statistical significance 
for both the — )■ J/4)D+ and — )■ J/iI'D*"'" signals estimated in this way is in excess of 
9 standard deviations. 

The low Q- value for the B+ — )• J/i[) D+ decay mode allows the B+ mass to be precisely 
measured. This makes use of the mass value, evaluated in Sect. |3} taking correctly into 
account the correlations between the measurements. The calibration of the momentum 
scale for the dataset used here is detailed in Refs. 27 33 . It is based upon large calibration 
samples of B+ — )■ J/i1jK+ and J/i|) — )■ |J."'"|J.~ decays and leads to an accuracy in the 
momentum scale of 3 x 10^^. This translates into an uncertainty of 0.30MeV/c^ on the 
B+ meson mass. A further uncertainty of 0.11 MeV/c^ arises from the knowledge of the 



detector material distribution [27,28,33,34 and the signal modelling. The uncertainty on 
the mass results in a 0.16 MeV/? uncertainty on the B^ meson mass. Adding these in 
quadrature gives 



6276.28 ± 1.44 (stat) ± 0.36 (syst) MeV/c^ 



The uncertainty on the D^^ meson mass and on the momentum scale largely cancels in the 
mass difference 



mr 



m 



4307.97 ± 1.44 (stat) ± 0.20 (syst) MeV/c^ 



5 Normalization to the B+ J/i|)7r+ decay mode 

A large sample of B+ — J/ipTr^ decays serves as a normalization channel to measure 
the ratio of branching fractions for the B+ — )■ J/\|; and B+ — )■ J/ij) 7T+ modes. Selection 
of B+ — > J/i[)7r+ events is performed in a manner similar to that described in Sect. [3] 
for the signal channel. To further reduce the combinatorial background, the transverse 
momentum of the pion for the B+ — t- J/ip 71+ mode is required to be in excess of 1 GeV/ c. 
The mass distribution of the selected B+ — )■ J/i|) 7t+ candidates is shown in Fig. |3} 

To determine the yield, an extended unbinned maximum likelihood fit to the mass 
distribution is performed. The signal is modelled by a double-sided Crystal Ball function 
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6.2 6.3 6.4 6.5 

m(J/\l)7r+) [GeV/c2] 

Figure 3: Mass distribution for selected — J/4)7r+ candidates. The results of a fit to 
the model described in the text are superimposed (solid line) together with the background 
component (dotted line). 

and the background with an exponential function. The fit gives a yield of 3009 ± 79 events. 
As cross-checks, a modified Novosibirsk function and a Gaussian function for the signal 
component and a product of exponential and polynomial functions for the background are 
used. The difference is treated as systematic uncertainty. 

The ratio of the total efficiencies (including acceptance, reconstruction, selection and 
trigger) for the B+ — t- J/ij) and — > J/i|) tC^ modes is determined with simulated data to 
be 0.148 ±0.001, where the uncertainty is statistical only. As only events explicitly selected 
by the JAIj triggers are used, the ratio of the trigger efficiencies for the B+ — )■ J/tJ; D+ and 
B+ — J/i|)7r+ modes is close to unity. 
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Figure 4: Distributions of (a) xiti^t) 0^) Xipi^t) fo^' ~^ J/U'Tt^ events: background 
subtracted data (red points with error bars), and simulation (blue histogram). 



6 Systematic uncertainties 

Uncertainties on the ratio TZj)+^^+ related to differences between the data and simulation 
efficiency for the selection requirements are studied using the abundant B+ — t- J/i|) 7t+ chan- 
nel. As an example, Fig. [4] compares the distributions of Xmi^t) Xipi^t) data 
and simulated — )■ J/ij)7t+ events. It can be seen that the agreement between data and 
simulation is good. In addition, a large sample of selected B+ — t- J/ip {K~^K~)^ K"*" events 
has been used to quantify differences between data and simulation. Based on the deviation, 
a systematic uncertainty of 1% is assigned. 

The agreement of the absolute trigger efficiency between data and simulation has 
been validated to a precision of 4% using B"*" — )■ J/i|; (K+K")^ K"*" events. A further 
cancellation of uncertainties occurs in the ratio of branching fractions resulting in a 
systematic uncertainty of 1.1%. 

The systematic uncertainties related to the fit model, in particular to the signal shape, 
mass and resolution for the B+ — )■ J/t|)D^ mode and the fit interval have been discussed in 
Sects. |4]and[5j The main part comes from the normalization channel B+ — )■ J/ipTT^. 

Other systematic uncertainties arise from differences in the efficiency of charged particle 
reconstruction between data and simulation. The largest of these arises from the knowledge 
of the hadronic interaction probability in the detector, which has an uncertainty of 2% per 
track p55]. A further uncertainty related to the reconstruction of two additional kaons 
in the B+ — )■ J/4)D+ mode with respect to the B+ — )■ J/4)7t"'" mode is estimated to be 
2 X 0.6% j36]. Further uncertainties are related to the track quality selection requirements 
Xtr < 4 and Pfake < 0.5. These are estimated from a comparison of data and simulation in 
the B^ — 7- J/4) 71"*" decay mode to be 0.4% per final state track. 

The uncertainty associated with the kaon identification criteria is studied using the com- 
bined B+ — )■ J/i[)D+ and B^ — )■ J/4)D*+ signals. The efficiency to identify a kaon pair with 
a selection on A^/" In C has been compared for data and simulation for various selection re- 
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quirements. The comparison shows a (—1.8 ±2.9)% difference between data and simulation 
in the efficiency to identify a kaon pair with 2 < min A^/" log L. This estimate has been con- 
firmed using a kinematically similar sample of reconstructed — t- J/i|) (K"'"K^)^ events. 
An uncertainty of 3% is assigned. 

The limited knowledge of the lifetime leads to an additional systematic uncertainty 
due to the different decay time acceptance between the B+ — )■ J/i[) D+ and B^ — )■ J/ij) 7r+ de- 
cay modes. To estimate this effect, the decay time distributions for simulated events are 
reweighted to change the B^ lifetime by one standard deviation from the known value |l], 
as well as the value recently measured by the CDF collaboration j2], and the efficiencies 
are recomputed. An uncertainty of 1% is assigned. 

Possible uncertainties related to the stability of the data taking conditions are tested 
by studying the ratio of the yields of B+ — )■ J/i|;K+7t+7r^ and B"*" — J/ijjK"'" decays for 
different data taking periods and dipole magnet polarities. This results in a further 2.5% 
uncertainty. 

The largest systematic uncertainty is due to the knowledge of the branching fraction 
of the D+ — )• (K~K"'")^j^ t& decay, with a kaon pair mass within ±20MeV/c^ of the known 
4) meson mass. The value of (2.24 ± 0.11 ± 0.06)% from Ref. [37] is used in the analysis. 
The systematic uncertainties on 7^^+^^+ are summarized in Table [3| 



Table 3: Relative systematic uncertainties for the ratio of branching fractions of — )■ J/ij^D^ 
and B+ J/i|j7r+. 



Source Uncertainty 

Simulated efficiencies 1.0 

Trigger 1.1 

Fit model 1.8 

Track reconstruction 2 x 0.6 

Hadron interactions 2 x 2.0 

Track quality selection 2 x 0.4 

Kaon identification 3.0 

B+ lifetime 1.0 

Stability for various data taking conditions 2.5 

E (d+ ^ (K-K+)^ 7t+ j 5.6 

Total 8.4 



The ratio '7^d*+/d+ estimated as 

^d:+/d+ = -jr- (2) 

where the ratio of yields is given in Table [2} The uncertainty associated with the 
assumption that the efficiencies for the B+ — )• J/4)D+ and B+ — )■ J/ipD*"*" modes are 
equal, is evaluated by studying the dependence of the relative yields for these modes for 
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loose (or no) requirements on the Xjp(B+), xiti^t) ^'■^d cr(B+) variables. For this selection 
the measured ratio of B+ — )■ J/i|;D*+ to B+ — )■ J/4'D+ events changes to 2.27 ± 0.59. An 
uncertainty of 4% is assigned to the '7^d*+/d+ ratio. 

The uncertainty on the fraction of the A±± amplitude, f±±, has been studied with 
different fit models for the parameterization of the combinatorial background, as well 
as different mass resolution models. This is negligible in comparison to the statistical 
uncertainty. 



7 Results and summary 



The decays B+ — t- J/i|)D+ and B+ — )■ J/il^D*"*" have been observed for the first time with 
statistical significances in excess of 9 standard deviations. The ratio of branching fractions 
for B+ — )■ J/t|>D+ and B+ — )■ J/i|)7t+ is calculated as 



i3(B+^ J/i|;7t+) B^^ " e 



B+^J/il)D^ 



A^(B+ ^ J/i|;7r+) 



(3) 



37 



with the mass of the kaon pair 



where the value of 8^^+ = B (^D+ (K-K+)^ tt^ 

within ±20MeV/c^ of the known value of the ^ mass is used, together with the ratio of 
efficiencies, and the signal yields given in Sects. |4] and |5j This results in 



g(B+^J/^|;D+) 
B{B+ J/i|jTr+) 



2.96 ± 0.67 (stat) ± 0.25 (syst). 



The value obtained is in agreement with the naive expectations given in Eq. (la)) and the 
values from Ref. [s] but larger than predictions from Refs. [7 10 



The ratio of branching fractions for the B+ — )• J/ip D*"*" and B+ — )■ J/i|) decays is 
measured to be 



g(B+^ J/i1;D: 
B{Bt^m^i 



2.36 ± 0.56 (stat) ± 0.10 (syst). 



This result is in agreement with the naive factorization hypothesis (Eq. (|lb[)) and with 
the predictions of Refs. [7|[8] . 

The fraction of the A±± amplitude in the B+ — )■ J/ijtD*"'" decay is measured to be 



r±± (B+^ J/x|;D; 
Etot (B+-^ J/i|;D* 



(52 ± 20)%, 



in agreement with a simple estimate of |, the measurements 



38 



39 



and factorization 



predictions |40j fo r B° — )• D* D*+ decays, and expectations for B^"— >■ J/ip^+A/^ decays 



from Refs. 41 , 42 
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The mass of the meson and the mass difference between the B+ and D+ mesons 
are measured to be 

= 6276.28 ± 1.44 (stat) ±0.36 (syst)MeV/c2, 
- = 4307.97 ± 1.44 (stat) ± 0.20 (syst) MeV/c^ 

The mass measurement is in good agreement with the previous result obtained by 
LHCb in the B^ — t- J/i|>7t'^ mode 28 and has smaUer systematic uncertainty. 
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